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Abstract. We describe all weighted Radon transforms on the plane for which the 
Chang approximate inversion formula is precise. Some subsequent results, including the 
r**- ! Cormack type inversion for these transforms, are also given. 



1. Introduction 

We consider the weighted ray transformation P\y defined by the formula 

43 ■ Pwf(s, e)= wise 1 - + te, 9)f( s e ± + te)dt, 

2 : R (1) 

s£t, 9 = (O u 9 2 ) e S 1 , 9 ± = (-6» 2 , 6»i), 

(N ■ i 
>• ■ where W = W(x, 9) is the weight, / = f{x) is a test function, x G R, 9 G § . Up to change 

of variables, Pw is known also as weighted Radon transform on the plane. 

We recall that in definition (1) the product RxS 1 is interpreted as the set of all oriented 

straight lines in R 2 . If 7 = (s, 9) G R x S 1 , then 7 = {x G R 2 : x = s9 ± + t9, t E R} 

(modulo orientation) and 9 gives the orientation of 7. 

We assume that 



W is complex — valued, 

W G C(R 2 x S 1 ) n L°°{R 2 n S 1 ), 

1 /• ( 2 ) 

io (x) d = — / W(x, 9)d9 ^0, x G R 2 , 
2n J 

where (i^ is the standard element of arc length on S 1 . 

If W = 1, then Pvk is known as the classical ray (or Radon) transform on the plane. 

If 

W(x,9) = exp(-Da(x,9)), 

+00 

r (3) 

Da{x,9) = / a(x + t9)dt, 



where a is a complex-valued sufficiently regular function on R 2 with sufficient decay at 
infinity, then Pw is known as the attenuated ray (or Radon) transform. 

The classical Radon transform arises, in particular, in the X-ray transmission tomog- 
raphy. The attenuated Radon transform (at least, with a > 0) arises, in particular, in the 
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single photon emission computed tomography (SPECT). Some other weights W also arise 
in applications. For more information in this connection see, for example, [Na], [K]. 

Precise and simultaneously explicite inversion formulas for the classical and attenuated 
Radon transforms were given for the first time in [R] and [No] , respectively. For some other 
weights W precise and simultaneously explicite inversion formulas were given in [BS], [G]. 

On the other hand, the following Chang approximate inversion formula for P\y, where 
W is given by (3) with a > 0, is used for a long time, see [Ch], [M], [K]: 



f apP r(x) = J h'(x9 ± ,9)d9, h'(s,9) = ^h(s,0), 

h(s,6) = -p.v. [ Pwf ^^ dt, s£l, 9 G S 1 , xe 1R 2 , 

7T J S ~t 



(4) 



where wq is defined in (2). It is known that (4) is efficient as the first approximation in 
SPECT reconstructions and, in particular, is sufficiently stable to the strong Poisson noise 
of SPECT data. The results of the present note consist of the following: 

(1) In Theorem, under assumptions (2), we describe all weights W for which the Chang 
approximate inversion formula (4) is precise, that is f appr = / on M 2 ; 

(2) For Pw with W of Theorem we give also the Cormack type inversion (see Remark A) 
and inversion from limited angle data (see Remark B). 

These results are presented in detail in the next section. In addition, we give also an 
explanation of efficiency of the Chang formula (4) as the first approximation in SPECT 
reconstructions (on the level of integral geometry) . 

2. Results 

Let 

d2\ 



C (R ) denote the space of continuous 
compactly supported functions on 

Let 



(5) 



L 00 ' "(R 2 ) = {/ : M a f G L°°(R 2 )} 
M a f(x) = (1 + \x\Yf(x), xeM 2 , a>0 



2 „ ( 6 ) 



Theorem. Let assumptions (2) hold and let f appr (x) be given by (4). Then 

f apP r = f (in the sense of distributios) on M 2 for all f G Co(M 2 ) (7) 
if and only if 

W(x,9) -w (x) =w (x) -W(x,-9), xeM 2 , 9 G S 1 . (8) 

(This result remains valid with C (R 2 ) replaced by L°°' fJ (M 2 ) for a > 1.) 
Theorem 1 is based on the following facts: 
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• Formula (4) coincides with the classical Radon inversion formula if W = 1 and, as a 
corollary, is precise if W = vjq. 

• Formula (4) is equivalent to the symmetrized formula 



fa PP Ax) = 47r J o(a;) J g\x9 ± ,9)d9, x e R 2 , 
< m 1 f Pwf(t,9) + P w f(-t,-9) ^ m clBv ei 



The following formula holds: 



(P w f(s,9) + P w f(s,-9)) = P Wsym f(s,9),(s,9) eRxE 1 , 
2 i (10) 

W sym (x, 9) = -(W(x, 9) + W(x, -9)), x e R 2 , 9 e S 1 . 



If 



g G C(R x S 1 ), suppq is compact, q(s,9) = q(—s, —9), 
g (s, 9) = -p.v. [ ^^-dt, (s, 9) G R x S 1 , 

7T J S-t 

R 

y g'(x9 ± , 9)d9 = as a distribution of 



(11) 



x G R 2 , 



then g = on R x S 1 . 

The statement that, under the assumptions of Theorem, property (8) implies (7) can 
be also deduced from considerations developed in [K]. 

Using that W sym = w under condition (8), we obtain also the following 

Remarks. Let conditions (2), (8) be fulfiled. Let / G C (R 2 ). Then: 

(A) Pwf on Q(D) uniquely determines / (or more precisely wof) on R 2 \D via (10) and 
the Cornack inversion from Pw sym f on ^(D), where D is a compact in R 2 , Q(D) 
denotes the set of all straight lines in R 2 which do not intersect D; 

(B) Pwf on R x (S + U S-) uniquely determines / on R 2 via (10) and standard inversion 
from the limited angle data Pw sym f on R x S + , where S+ is an arbitrary nonempty 
open connected subset of S 1 , S- = —S + . 

For the case when W is given by (3) under the additional conditions that a > 
and suppa C D, where D is some known bounded domain which is not too big, and for 
/ G C(R 2 ), / > 0, supp f C -D, the transform Pwf is relatively well approximated by 
p w appr f, where W appr (x, 9) = w (x) + (l/2)(W(x, 9) - W(x, -9)). In addition, this W appr 
already satisfies (8) . This explains the efficiency of (4) as the first aproximation in SPECT 
reconstructions (on the level of integral geometry) . 
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